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ABSTRACT 

Spinning planetesimals in a gaseous protoplanetary disk may experience a hydrody- 
namical force perpendicular to their relative velocities. We examine the effect this force 
has on the dynamics of these objects using analytical arguments based on a simple 
laminar disk model and numerical integrations of the equations of motion for individ¬ 
ual grains. We focus in particular on meter-sized boulders traditionally expected to 
spiral in to the central star in as little as 100 years from 1 A.U. We find that there are 
plausible scenarios in which this force extends the lifetime of these solids in the disk 
by a factor of several. More importantly the velocities induced by the Magnus force 
can prevent the formation of planetesimals via gravitational instability in the inner 
disk if the size of the dust particles is larger than of order 10 cm. We find that the 
fastest growing linear modes of the streaming instability may still grow despite the 
diffusive effect of the Magnus force, but it remains to be seen how the Magnus force 
will alter the non-linear evolution of these instabilities. 

Key words: Planetary systems - protoplanet ary disks, Planetary systems - mete¬ 
orites, meteors, meteoroids 


1 INTRODUCTION 

Planets are common, but the process by which they form 
may be complicated and is difficult to observe directly. 
Somehow the small dust grains present in the interstellar 
medium of galaxies must become > 1000 km-sized planets. 
Many orders of magnitude in this ‘size ladder’ are reasonably 
well-understood. Below ~ 10 cm, grains can grow by stick¬ 
ing together during collisions (see Blum &; Wurm 2008, for 
an extensive review of this collisional process), while grav¬ 
ity is strong enough for objects larger than a few km that 
they can grow by gravitationally-focused collisions. Between 
these regimes lies the meter barrier, where solids are both 
too small for their self-gravity to be important, and are likely 
moving too fast relative to the gas and hence each other to 
stick upon colliding. To make matters worse, meter-sized 
particles experience particularly effective drag forces, caus¬ 
ing them to spiral into the star at the center of the disk in 
an astronomically minuscule time of order 100 years from 
1 A.U. It remains a mystery how exactly the solids which 
eventually become planets overcome this barrier. 

Among the most influential ideas is direct gravita¬ 
tional collapse through the collective gravity of many grains. 
Goldreich & Ward (1973) presented an early version of this 
idea, in which dust would settle out of the gas disk into 
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a dynamically cold midplane where the dust disk would 
be subject to Toomre instability leading to direct col¬ 
lapse to ~ 10 km planetesimals. Although this scenario 
is unlikely owing to turbulence in the disk and Kelvin- 
Helmholz-like instabilities which increase the grain veloc¬ 
ity dispersion (Weidenschilling 1980), a variety of physi¬ 
cal mechanisms have been proposed to aid the grains in 
their gravitational collapse, including pressure traps as¬ 
sociated with magneto rotational instability-induced tur¬ 
bulence (Fromang & Nelson 2005) and the streaming in¬ 
stability in which the drag on many individual parti¬ 
cles has a large back-reaction on the gas when the den¬ 
sity in solids is of order the density of the gas (e.g. 
Goodman & Pindor 2000; Youdin & Goodman 2005). A re¬ 
cent review by Johansen et al. (2014) details these pro¬ 
cesses and other aspects of planetesimal formation. 

In this paper we explore the effects of another poten¬ 
tially important piece of physics. First described by Newton 
(1671), the Magnus force is best-known for its effects on ter¬ 
restrial sports (Mehta 1985) and aeronautics (Seifert 2012). 
Spinning objects moving with respect to a background fluid 
create an asymmetric wake in the fluid around them, which 
in turn alters the forces felt by the object. In addition to the 
drag force which opposes the motion of individual grains rel¬ 
ative to the background gas, a component of the fluid’s force 
on the object is perpendicular to the relative motion. When 
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the circulation is caused primarily by the spin of the object, 
the force will also be perpendicular to the spin axis. 

The addition of a force beyond the usually-considered 
drag and gravitational forces has the potential to change the 
mechanics of many of the scenarios commonly considered in 
planet formation, particularly for small objects where in¬ 
teraction with the gas is already known to be important. 
We can immediately guess that the force will tend to be 
diffusive, in that it will apply forces to individual particles 
isotropically, assuming that whatever process gives rise to 
the spin does not give rise to a preferred spin axis. All else 
being equal, we therefore expect that this effect will make 
planetesimal formation more difficult. 

We will argue that there are plausible physical mecha¬ 
nisms that provide meter-sized solids with enough spin for 
the Magnus force to play an appreciable role in the dynamics 
of these objects. We will also argue based on our simulations 
of the dynamics of individual grains that in regimes where 
the Magnus force acts at all, disks of solid material are un¬ 
likely to be dynamically cold enough to collapse gravitation¬ 
ally in the inner disk unless the grains undergoing collapse 
are smaller than about 10 cm. 

In section 2, we collect the basic ingredients of a laminar 
disk model and its interactions with individual particles. We 
detail our treatment of the Magnus force in section ??, then 
explore the conditions under which it is relevant to the dy¬ 
namics in section 4. Focusing on one such regime, we carry 
out numerical integrations of the equations of motion includ¬ 
ing the Magnus effect in section 5. In section 6 we analyze 
the results in a co-rot at ing reference frame, and show that 
the dynamics we observe in the integrations can be largely 
explained as the particles following force-free trajectories in 
this frame. We use this fact to map out the effects of the 
Magnus force in several planetesimal formation scenarios as 
a function of heliocentric location in the disk, particle size, 
and other relevant parameters. We summarize in section 7. 


2 CLASSICAL SCENARIO 

In this section we review the key issues regarding the flow 
of planetesimals in a proto-planetary disk. For simplicity 
we adopt a simple disk model following Chiang & Youdin 


( 2010 ), 

wherein 
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respectively the gas density, gas scale height, gas mean free 
path, fractional velocity suppression of the gas by pressure 
support, gas sound speed, and gas temperature. We also 
define the Keplerian velocity vk = tQk = yjGM*/r, for a 
central star of mass M*. The cylindrical distance from the 
star is given in astronomical units as tau, while ks and 
rriH represent Boltzmann’s constant and the atomic mass of 


hydrogen. F is a factor by which the disk’s overall density 
may be scaled. 

From the numbers above, we have almost all the infor¬ 
mation necessary to calculate the grain Reynolds number 


Re — 


4 SV re l 


Ac, 
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where s the radius of the grain. To find v re i , the relative 
velocity between the grain and the gas, we need to take 
into account the drag force felt by the grain. In the limit of 
ineffective drag, the magnitude of this velocity v re i ~ t/vk , 
i.e. the planetesimal will move at nearly Keplerian velocity, 
whereas the gas will be slightly slower owing to its pressure 
support. Small grains that are well-coupled to the gas will 
have substantially smaller velocities, however. 

The drag force experienced by individual grains oper¬ 
ates in one of four regimes depending on the Reynolds num¬ 
ber and the mean free path of gas molecules relative to the 
grain size. We take the dimensionless drag coefficient Cd, 
defined as the ratio of the drag force FT to (l/2)p g 7vs 2 vl eh 
to be 


C D = 
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Stepinski & Valageas (1996). See 

also 


following 

Weidenschilling (1977); Garaud et al. (2004); Loth (2008). 
In the high-Re regime where inertial forces dominate 
viscous forces, the drag arises when the oncoming gas is 
deflected and slowed by the particle. In the high-A regime 
where the fluid equations break down, the drag arises via 
a slight difference in the momentum flux on either side of 
the particle, in part owing to the thermal motion of the 
gas particles. In between, viscous forces become important 
increasing the drag relative to the inertial regime. 

With the drag force in hand, we may now define the 
stopping time t s = mv re i/Fd and its dimensionless version 
t s = ts^K- The velocity of the particle relative to Keplerian 
may be approximated following Chiang & Youdin (2010), 


and 


v r « 2pQ K r 

1 + r s 2 
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See also Weidenschilling (1977); Youdin (2010). The particle 
has a tangential velocity relative to the gas of — (1 — 
rj)vK — tjvk + — vk), yielding a total relative velocity 

magnitude 


_ o r sV 4 + r| , . 

Vrei ~ rjll K r —; 5 — (11) 

1 + T£ 

We can now plug in the definition of r s , which depends in 
a non-trivial way on v re i , and solve numerically for v re i , 
which closes the whole system of equations, letting us find 
Re, Fd, and t s . Figure 1 shows the resulting dimensionless 
stopping time. Small particles close to the center of the disk 
have t s C 1 , meaning they are strongly coupled to the gas, 
with velocities closer to the gas velocity than the Keplerian 
velocity. Large rocks far from the star are weakly coupled, 
and move at very nearly the Keplerian velocity. 
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Figure 1. The log of the non-dimensional stopping time of grains 
of various sizes at various locations in the disk. We show the 
contour where r s = 1, i.e. the dividing line between the grains 
being well-coupled to the gas via drag, and the grains largely 
ignoring the gas. 



log 10 (r) (AU) 


Figure 2. The log of the inflow timescale r/v r given in years. We 
highlight 100 yr, 1000 yr, and 1 Myr. In the classical scenario, few 
if any particles can survive the full multi-Myr lifetime of the disk, 
and objects of order 1 meter in size can quickly spiral in through 
the disk to the star. 


Since we know the relative velocities, we can evaluate 
the drag regime in which these particles live by checking the 
conditions in equation 8, and their inward velocity v r with 
equation 9. We can therefore estimate the timescale on which 
the particles will lose a substantial fraction of their orbital 
radius as ti n — r/v r , shown in figure 2. Particles of order 1 
meter in size, spiral in to the star in of order 100 years from 
1 A.U. This is the classical meter barrier problem. 


3 THE MAGNUS FORCE 

The lift force experienced by spinning objects has been stud¬ 
ied analytically, numerically and experimentally over the 
past century. A typical case involves a sphere or cylinder 
spinning about its axis of symmetry suspended in a wind 
tunnel, with the background velocity perpendicular to the 
spin axis. The resultant forces on the cylinder can be mea¬ 
sured (e.g. by measuring the tension in the apparatus sus¬ 
pending the cylinder) or computed (by integrating the ap¬ 
propriate components of the stress tensor along the surface 
of the object). 

Just as with the drag force, the lift force Fl may then 
be expressed as a dimensionless coefficient, 

Cl = a (12) 

(1/2)7 TPgVlelS 2 

which should in principle only depend on the dimensionless 
parameters of the flow, Re, S = suj/v re i , and A/s. Unlike the 
drag force, which always points towards —v r eh the lift force 
refers to any force on the object perpendicular to D re z, so it 
is in principle a 2D vector quantity. Though some authors 
do indeed investigate it this way (e.g. Poon et al. 2013), for 
simplicity we take Cl to be a scalar quantity with ft m 
uj x Vrei/\w x v re i\, so that positive (negative) values of Cl 
refer to lift forces aligned (anti-aligned) with uj x v re i. When 
the spin and velocity axes are aligned, we take Cl = 0. By 
symmetry we take the lift force to be entirely parallel or 
anti-parallel to n, though clearly this will not always be the 
case. 

Among the best-known results in the long history of 
the lift force is that of Rubinow & Keller (1961), who an¬ 
alytically computed the flow around a spinning sphere in 
the limit that Re <C 1. They derived a value for the lift 
coefficient 


^Rubinow^^^! ( 13 ) 

to leading order in Re. This is a reasonable approximation 
to Cl when Re< 1, so long as A/s < 4/9. This turns out 
to be a rare circumstance in our fiducial disk model, in that 
once Re< 1, typically A/s >4/9. When A/s >4/9, we as¬ 
sume that Cl =0, since the fluid no longer ‘cares’ about 
the physical size of the object. This means that Cl will 
rapidly drop from arbitrarily large values to zero for Re ~ 1 
and large S. The nature of this transition is uncertain and 
not constrained by data. 

At higher Re, no analytic results are available. Individ¬ 
ual studies will typically run a series of laboratory or nu¬ 
merical experiments over a small dynamic range of Re and 
S. Many report fitting formulae for their results, though 
of course experimental error and the choice of functional 
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form can make these unreliable even over the range of pa¬ 
rameter space probed by the experiments. For instance, 
Oesterl‘e & Bui Dinh (1998) and You et al. (2003) cover a 
similar range in Re and <S, but their formulae disagree by as 
much as 50% for Re « 20 and S « 1. 

An extensive review by Loth (2008) has assembled a 
large quantity of historical data on the lift owing to both 
particle spin and background shear in the fluid flow. He pro¬ 
poses the following global fit for Cl as a function of Re and 
S for low to moderate Re< 1000. 

Cp th = 25 [1 - (0.675 + 0.15( 

1 + tanh(0.56(5 — 1)))) tanh(0.18JJe 1 ^ 2 )] (14) 

This fit encompasses data from Legendre & Magnaudet 
(1998), Bagchi & Balachandar (2002), Tri et al. (1990), and 
Tsuji et al. (1985), and approaches the Rubinow & Keller 
(1961) solution in the limit Re —>• 0. 

At higher Re , the picture becomes murkier owing to 
the onset of turbulence in the wake of the sphere. For 
Re >40, 000, Loth (2008) describes two regimes, both with 
values of Cl wildly divergent from the formula quoted above 
when S < 1. In the subcritical regime, where the surface of 
the object is smooth enough that the boundary layer re¬ 
mains laminar as it separates from the object, many authors 
have noted that Cl changes sign (Macoll 1928; Davies 1949; 
Tani 1950) at low S. Kim et al. (2014) have recently de¬ 
veloped a detailed model that predicts where and how this 
‘inverse’ Magnus effect occurs based on the results of wind 
tunnel measurements for 0.6 x 10 5 < Re < 1.8 x 10 5 and 
0 < S < 2 . We expect that this regime is irrelevant for solids 
in protoplanetary disks owing to the irregular shapes in¬ 
volved, analogous to the dimples on golf balls (Davies 1949; 
Bearman & Harvey 1976). 

When the particles are rough, or Re is sufficiently 
large, turbulent rather than laminar separation occurs. In 
this ‘super-critical’ regime, observational data suggest that 
Cl can become quite large again, nearly approaching the 
Rubinow & Keller (1961) limit. Fitting to a series of mea¬ 
surements collected using baseballs, Sawicki et al. (2003) 
suggest 

^Sawicki = £ min ( L5) 0 6 + 0.09/5) (15) 

which agrees well with data from Nathan et al. (2006), 
Briggs (1959), and Watts & Ferrer (1987) for Re > 10 5 . Un¬ 
fortunately, it is not obvious how one would connect the low 
Cl at moderate Re to these high values of Cl at high Re. 
Moreover, the Sawicki et al. (2003) fit over-predicts Cl for 
high S where the separation between sub- and super-critical 
regimes is expected to disappear. For the sake of continuity 
with the moderate-Re regime, we therefore adopt 

Tanaka min ( 0 . 5j 0 . 5 / 5 ) (16) 

from Tanaka et al. (1990) in the high-Re regime. This is a 
conservative choice, and it is possible that in this regime we 
may be underestimating the Magnus force by about a factor 
of two. 

We thereby arrive at the following formula for the lift 
coefficient, which we expect to be a reasonable approxima¬ 
tion for objects in protoplanet ary disks over at least 5 orders 


of magnitude in Re and S < 10, 

c = I |o> X Vrei I max(C'p th , C£ anaka ) for A/s < 4/9 
1^0 for A/s > 4/9 

(17) 

with the caveat that the model of Kim et al. (2014) and the 
fit from Sawicki et al. (2003) may be relevant for S < 1 and 
Re > 10 5 . 

The ratio of the lift force to the drag force, assuming 
that the lift is perfectly aligned with n, is simply Cl/Cd- 
This quantity should give us an idea of when and where the 
lift force is important. The three limiting cases (neglecting 
the orientation factor) are 

{ Y 2 <S Re Re 1 and A/s < 4/9 
1.15 5 <1 and Re > 1 (18) 

1.1 S 1 and Re 1 

Excluding the possibility of extremely large dimensionless 
spins, the maximum values of Cl/Cd are likely to be found 
in the high-Re regime, which tends to correspond to the 
inner regions of the disk and to larger grains. 


4 THE SPIN RATE 
4.1 Empirical values of uj 

We now turn to the question of what value to use 
for uj. There is a substantial quantity of data available 
(Warner et al. 2009), derived largely from photometric light 
curves, regarding the spin states of modern asteroids with 
sizes s > 10 meters, though of course these asteroids have 
been subject to several Gyr of evolution under a different 
set of environmental conditions than existed in the early so¬ 
lar system. These data show (see Figure 3) that asteroids 
larger than 100 m have a minimum period of a few hours, 
consistent with the idea that they are largely held together 
gravitationally, in which case the minimum period is roughly 
a/37t/p s G, where p s is the mean density of the body. Below 
100 meters, there are plentiful rapidly-rotating asteroids, in¬ 
cluding some with periods of a minute or less. With this in 
mind we adopt the following functional form 

{ ^0 ( s ,, , 1 S ^rubble / in \ 

\ rubble/ (19) 

^0 ^ ^rubble 

Above a critical size s ru bbie we take particles to be rubble 
piles held together by gravity with a constant rotation period 
uo of order the free fall time. Smaller particles are taken to 
have a power law index —/3 with smaller grains rotating 
faster. 

The black and red dashed lines in figure 3 show exam¬ 
ples with (cco, Srubbie, yd) = (3 x 10 -4 seconds -1 ,1 km, 1.0), 
the shallow case, and (3 x 10 -4 seconds -1 ,0.3 km, 1.5), the 
steep case, respectively. These lines are plausible, though 
simplified, representations of the data. We are mostly con¬ 
cerned with much smaller-sized objects, so even if we be¬ 
lieved the size-frequency distribution of solids in the solar 
system to be the same today as it was in the presence of a 
gaseous disk, these relations would still be huge extrapola¬ 
tions owing to the difficulty of observing small solar system 
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Figure 3. The observed angular velocity-size relation of modern 
asteroids from Warner et al. (2009, February 2014 version). The 
logarithmically-spaced contours show the density of the highest- 
quality (U > 2) data points in the catalog where both parameters 
have measured values, rather than just upper or lower limits. The 
black and red dashed lines show equation 19 evaluated with the 
same wo = 3x 10 -4 seconds -1 , but with two different power law 
slopes (—1 and —1.5) and two different cutoffs (1 km, 0.3 km). 


bodies. There is other evidence that small meteor-sized ob¬ 
jects may rotate rapidly, namely the appreciable fraction of 
non-linear meteor trails (Beech 1988). 

It is worth pointing out that if f3 > 1, eventually the 
grains will reach their theoretical maximum spin rate of 
ujs 1 km s-\ the sound speed of solid rock. This will 
happen for s < s r ubbie(^oSrubbie/l km s -1 ) 1 /^ -1 ^. In subse¬ 
quent figures, we cap ujs at 0.01 km s _1 , which should be 
substantially slower than the maximum spin rate, though 
still perhaps implausibly large. 

For these values of uj as a function of the object’s size, 
we can estimate the ratio of the Magnus force to the drag 
force as a function of quantities we know how to calculate 
in the classical scenario. 

Clearly this ratio depends quite sensitively on what we 
assume about the size-frequency distribution - the Magnus 
force may be small but not negligible in one case, or over¬ 
whelmingly large in the other. Figure 4 shows the ratio for 
the shallow cj(s) case, while 5 shows the steeper case. When 
uj oc s -1 , the features in this graph come mostly from v re i - 
when the relative velocity of the particle and the gas is min¬ 
imized, Fl/Fd is maximized. This effect can be completely 
swamped by any change in the o;(s) distribution, as in figure 
5, where rapidly rotating well-coupled (r s 1) grains show 
the most prominent influence of the Magnus force. 


4.2 Theoretical estimates for uj 

Given the uncertainty in using the spin rate distribution 
^empirical > it is also worth considering the physical mecha¬ 
nisms responsible for the spin. This has its own set of uncer¬ 
tainties, including unknowns in the basic physics, and the 
sensitivity some of these physical mechanisms have to un¬ 
known quantities, such as the small-scale shapes of the solid 



logio(r) (AU) 


Figure 4. The log of the ratio of the Magnus force to the drag 
force, with contours at —1 and —0.5. The spin rate decreases 
as s -1 , consistent with data for modern asteroids. The Magnus 
force under this scenario is not particularly effective for meter¬ 
sized objects, but may be promising for small grains close to the 
star. The white region at large radii corresponds to Re < 1, where 
we believe the Magnus force is unlikely to play a role. 


bodies, or the instantaneous distribution of particle spins, 
sizes, and orbits. Nonetheless we discuss a few potentially 
important ingredients. 


4.2.1 YORP 


A great deal of effort has gone into explaining the spin of 
modern asteroids with the YORP effect (Rubincam 2000; 
Vokrouhlicky & Capek 2002), caused by the anisotropic re¬ 
radiation of incident sunlight. Theoretically this torque 
scales roughly as 


Ty 


BGi 

4 a 2 Vl — e 2 p s s 2 


Cyl 


( 20 ) 


following (Rossi et al. 2009). The orbital semi-major axis 
and eccentricity, a and e appear in the denominator, mak¬ 
ing the force more effective near the center of the disk for 
particles on eccentric orbits. The solar radiation constant G i 
is taken to be 10 22 g cm s -2 , while B « 2/3 is the Lamber¬ 
tian scattering coefficient. The coefficient Cy is a number 
variously measured for solar system asteroids to be between 
-0.025 and 0.025 (Rozitis & Green 2013), while the moment 
of inertia I « (tt/2 )p 3 s 5 . 

We therefore expect the timescale on which a particle’s 
angular velocity may be changed by YORP is 


£yorp ~ ujI/\Ty\ ~ 10 6 cJorAi/ps,oSo seconds (21) 


where we have adopted Cy — 0.01, and normalized uj, p s 
and s to unity in cgs. Modern asteroids exhibit Cy of both 
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i°gioM < AU ) 

Figure 5. The log of the ratio of the Magnus force to the drag 
force, with contours at —1 and —0.5. The spin rate decreases as 
s -1 - 5 , again consistent with data for modern asteroids. This case 
shows similar behavior to the shallower u(s) distribution, though 
here the steeply increasing spin rate for smaller grains means that 
grains smaller than 10 cm are strongly affected by the Magnus 
force, provided they have a high enough Re for such fluid effects 
to be important. 


signs, meaning YORP can spin up an asteroid (which if not 
countered by some other torque may lead to its disruption 
(Jewitt et al. 2014; Jacobson et al. 2014)) or spin it down to 
near-zero u, perhaps explaining the population of extreme 
slow rotators. 

YORP has the disadvantage of requiring direct irradi¬ 
ation by the star, which may be impossible for much of the 
disk. A gaseous disk with homogeneously-mixed 0.1 pm dust 
can be quite optically thick (r ~ 10 5 ) even in the vertical 
direction at optical wavelengths (Chiang & Goldreich 1997). 
In order for YORP to be effective, these small grains would 
need to be depleted in favor of much larger grains, since the 
opacity in the geometric limit scales as the surface area to 
volume ratio of an individual grain. 


4-2.2 Hydrodynamic friction 

A torque unlikely to be important in today’s solar system 
due to the absence of gas, but likely quite important during 
the disk’s lifetime, arises from drag along the surface of the 
body. It is not precisely clear how to estimate this torque, 
as it may well depend on the shape and small-scale features 
of the object, not to mention the detailed flow of the gas. 


As a rough estimate, we take 


X 

— \u X V re [ | 

(22) 

Ftop 

= eF D (v = \v re lX + us\) 

(23) 

Fbot 

= eF D (v = \v re lX - tdS |) sign {us - Vreix) 

(24) 

Td 

— — s\Ftop + Fbot |/2 

(25) 


Here x is just the fraction of the velocity that is not aligned 
with the spin axis, and should therefore be included in the 
relative velocity between the object and the surrounding gas. 
F t0 p and Fbot are the tangential forces exerted on the skin of 
the object on the side rotating in to the oncoming relative 
velocity and the side heading away from that velocity re¬ 
spectively. These are estimated to be some fraction e of the 
overall drag force the object as a whole would feel if it were 
moving at the surface velocity relative to the gas. The sign 
change when su = v re i comes about because in the limit of 
fast rotation, both forces will oppose that rotation, whereas 
for slow rotation the forces should nearly cancel each other 
out, until at zero rotation this drag torque is zero. 

In the limit where the particle is spinning rapidly, i.e. 
S 1, we find that Tb = —eCDp g 7rs 5 u 2 , where of course 
Cd depends on the velocity us and location in the disk via 
the Reynolds number. We can compare this to an expression 
from Loth (2008) for the torque in the limit where v re i —>■ 0 
and the background shear is zero, 

T D ,Loth = -Kl6ps 6 uj 2 ReF (l + . (26) 

We use Re^ to denote the Reynolds number where the char¬ 
acteristic velocity is us instead of v re i. This agrees with 
the expression for T given by Rubinow & Keller (1961) for 
Reu <C 1. In this regime, we know that Cd — 24/ Re, so the 
two expressions agree for e = 2/3. This gives us some confi¬ 
dence about our expression for Td, although it still remains 
an extrapolation to apply it when S < 1. 

In the high-spin regime, the spindown time is 

^spindown,fast = ^1 /\Td\ ^ ™ (27) 

pg ZCCd 

In other words the particle will spin for some large multiple 
of its rotational period, of order or a bit larger than the den¬ 
sity contrast. In the opposite limit S <C 1, the drag torque 
is a little more complicated according to equation 25, 

Yd, slow ^ cCd'F'X PgS VrelU ^1 Y “ I v=v rel ^ (28) 

Note that in general the logarithmic derivative 
will be between —1 and 0 depending on the 
drag regime. For convenience we will abbreviate 
C = (1 + (l/2)d\nCD/d\nv\ v=Vrel ), which will always 
be of order unity. In this limit the spindown time is 

^spindown,slow — ^spindown,fast X (29) 

i.e. substantially shorter than in the fast-spinning limit. This 
is somewhat counter-intuitive because in this limit the two 
drag torques are acting in opposite rotational directions in¬ 
stead of acting together to slow down the rotation. However, 
the drag in this regime is much more effective (all other 
parameters being equal), because the drag forces F top and 
Fbot care about v re i , which is large compared to the velocity 
which is being decreased, us. 
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Comparing equation 21 and equations 27 and 29 we see 
that each timescale has a different dependence on ca, with the 
YORP timescale increasing with increasing a;, but the spin 
down times decreasing or remaining constant with increasing 
uj. This means that if these two forces were to counteract 
each other, they would pick out an equilibrium u, which we 
can find by setting the timescales equal. It turns out that 
in this equilibrium, the slow-rotation limit is appropriate, in 
which case we obtain 


r At / 7 (js 1 + Ts)/ V 4 + 

6.65 2eC D s oX 2 C 


(30) 


Remarkably, if we focus on the high-Re and loosely coupled 
(r s 1) regime, and evaluate for a = 3 AU, appropriate 
for the modern asteroid belt, we obtain ujo ~ 2.6/ex 2 s o, 
numerically very similar to the shallow ^empirical case when 
X - e - 1. 

We have so far assumed that surface drag will mono- 
tonically spin down a particle, but this is not guaranteed. 
For particles with asymmetries derived from an empirical 
model, Capek (2014) has found that appreciable spin can be 
obtained by the action of hydrodynamic drag in the context 
of comet ary meteoroids. The median spin rate for particles 
spun up in this manner is 


S ~ 


MT 3 # 12 . 


(31) 


The deviation from S ~ const, arises from the assumed scal¬ 
ing of asymmetries in particle shape as a function of size. 
The simulations are also carried out in a low-Re regime 
where we expect the Magnus force to be unimportant, 
Nonetheless, this may provide an approximate lower limit 
on the spin rate. Moreover, although the median is rather 
small, the distribution is quite broad - log-normal with 0.5 
dex scatter. Particles in the simulations in question can even 
reach breakup spin rates. 


4-2.3 Collisions 

The gravitational torques exerted by the star (or planets) 
on planetesimals, commonly included in evolutionary mod¬ 
els for the asteroid belt, are probably unimportant given the 
long timescales over which they act and the smallness of the 
objects in which we are most interested. Similarly, accretion 
torques and shear of the disk across the size of the solid 
body should be small effects. Collisions, however, may be 
important or even dominant in determining the spin rate dis¬ 
tribution, but modeling them may require a self-consistent 
evolution of a spatially-dependent size-frequency distribu¬ 
tion. Nonetheless we can take some straightforward steps in 
that direction. 

We begin by adopting some idealizing assumptions 
which can be relaxed later. We assume that all particles 
are the same size with an isotropic velocity dispersion <j, 
local number density n, and cross-sectional area ns 2 . In 
this case the characteristic time between collisions for a 
particular particle is t c ~ (nans 2 ) -1 . If the particles are 
in a thin disk, we can take the velocity dispersion to be 
a = vkH/v = QkH, where H is the scale height of the par¬ 
ticles. The column density of particles (number of particles 
per unit area) is just N s = nH : so t c ~ (jV s 7rs 2 fIic) _1 . By 
assuming all particles to be the same size, we can estimate 


the column density of particles to be the mass surface den¬ 
sity of solids divided by the mass of an individual particle, 
i.e. N s = 3T, s /(4np s s 3 ). Adopting the mass surface density 
of solids 

= 33FZ re ir^j /2 g cm -2 (32) 

from Chiang & Youdin (2010), where Z re i is the metallicity 
relative to solar, we find the following typical time between 
collisions, 

t c ~ 0.04(FZ re 0 _ V 3/ ^s o p s ,o/^K. (33) 


In other words, particles can collide more frequently than an 
orbital time for small particles near the center of the disk. 

Perhaps more important than the orbital time are the 
other timescales on which we expect the spin may change. 
We can envision a scenario in which collisions at some ve¬ 
locity v c set a maximum spin rate of uj c ~ v c /s, which then 
decays on a spindown timescale. Particles can reach this 
maximum spin rate following a collision at a large impact 
parameter with an object of similar size, so long as a sub¬ 
stantial fraction of the impactor’s angular momentum in the 
target’s reference frame ends up as spin angular momentum 
of the target (Farinella et al. 1992). For smaller fragments 
of a collision, rapid rotation is also a common outcome (e.g. 
Fujiwara &; Tsukamoto 1981; Paolicchi et al. 1989). 

The typical spin rate for objects subject to such col¬ 
lisions will then either be uj c if t c <C t sp i n down , or a much 
smaller value set by some other process (e.g. the YORP- 
spindown balance posited above) when t c t sp indown- The 
ratio of these two timescales evaluated at uj = uj c is roughly 


tc cCd 3/14 

7 - ~ry V C,3^' AU 

^spindown ^rel 


(34) 


where we have normalized the collision velocity to 10 3 cm/s. 
Enhanced particle densities, as commonly present in direct 
gravitational collapse or streaming instability scenarios, will 
decrease this ratio in direct proportion to the overdensity. 
We therefore consider it quite plausible that collisions may 
become frequent enough over the course of the planet forma¬ 
tion process that even effective gas drag would be insufficient 
to slow the particles much below their maximum spin value 
(jj c - If we make the further approximation that v c ~ v re i , 
the dimensionless spin rate would approach S ~ 1. In this 
scenario, equation 18 tells us that the ratio of the lift to the 
Drag force would be of order unity for grains larger than the 
mean free path of gas particles in the disk. 


5 NUMERICAL INTEGRATIONS 

We have shown that solid particles in our fiducial protoplan¬ 
etary disk can plausibly have spins large enough for the lift 
force and the drag force to be comparable. The next step is 
to understand the effect this new force has on the dynamics 
of single grains. To that end, we carry out explicit numer¬ 
ical integrations of the equations of motion for individual 
particles traveling through the simple laminar disk model 
presented in section 2. 

In section 4, we saw that there are many plausible ways 
to assign a spin cj to particles in the disk, depending on their 
size, location, and most importantly the physical processes 
that are relevant in determining the spin. In order to carry 
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Figure 6. The semi-major axis as a function of time for 100 particles with randomly-oriented spins. Each panel shows a different value 
for the initial spin rate, while the color indicates how well the spin-axis is aligned with the orbital angular momentum. The black line 
shows the case with zero spin for reference. Successively faster and more anti-aligned spins both tend to increase the lifetime of a particle 
in the disk. For very rapid spins, the orientation starts to matter much less, and all particles survive much longer than in the reference 
case. 
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Figure 7. Various properties of the integrations for *Starget = 1. 
As in Figure 6, the colors represent the alignment of the spin with 
the orbital plane, and the no-spin case is shown as the black line. 
From top to bottom, the quantities are inclination, relative veloc¬ 
ity, Magnus force ratio, dimensionless stopping time, and the nor¬ 
malized azimuthal component of v re i . The Magnus force induces 
oscillations in all of these quantities, and substantial departures 
from the no-spin case. 


Figure 8. Various properties of the integrations for ^target = 1. 
Same as Figure 7, except that the quantities have been time- 
averaged in 1-year windows and each quantity is shown as a func¬ 
tion of semi-major axis instead of time. Particles therefore tra¬ 
verse this plot from right to left. This smooths out the oscillations 
seen in the previous figure, and accounts for the fact that at a 
given time the particles will be at different locations in the disk, 
and hence subject to different densities, gas velocities, etc. The 
solid black line again shows the zero-spin trajectory. 
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out integrations, we need to make an explicit choice. The 
simplest approach would be to pick a constant value of uo 
at the beginning of the simulation and leave it unchanged 
throughout. However, since several scenarios pick out a par¬ 
ticular dimensionless spin <S, rather than a dimensional spin 
cj, we set the spin rate such that S will exponentially ap¬ 
proach some constant <Starget, and we run several sets of 
simulations with different values of <S ta r g et • 

We also expect that the orientation of the spin axis 
Cj will evolve over time. Each time a particle collides, or 
is subject to enough time under hydrodynamical or radia¬ 
tive torques, the spin may be reoriented. However, we have 
found that keeping the spin axis constant in time makes 
interpreting the results of the simulations more straightfor¬ 
ward. This is because, as we shall see, the particle’s behavior 
is largely determined by lj z , the projection of the spin orien¬ 
tation on the angular momentum vector of the gaseous disk. 
One should therefore keep in mind that we do not expect 
any particular particle to follow these exact trajectories. In¬ 
stead, each particle behaves according to its instantaneous 
semi-major axis and spin orientation. 

We evolve each particle’s position, velocity, and spin in 
an inertial frame with Cartesian coordinates x, y, and z. The 
disk is centered at the origin with its angular momentum 
parallel to z. In this system, each particle is subject to the 
following 7 coupled ODE’s 


dx 

dt 

dv 

dt 

d\u\ 

dt 


v 



+ 


Fp + Fl_ 

rris m s 


(target -5)-LL 
t-orb 


(35) 

(36) 

(37) 


Here the particle mass m s = (47r/3 )p s s 3 with p s = 

1 g cm -3 . The particle’s spin is set to approach a tar¬ 
get dimensionless spin <Star g et on a local orbital timescale 
torb = &K 1 - For the sake of simplicity and clarity, the orien¬ 
tation of the particle’s spin and the particle’s size are taken 
to be constant throughout the simulation. 

The drag and Magnus forces, Fp and Fl are calculated 
through equations 8 and 17. Each depends on u re z, which is 
calculated instantaneously as the vector difference between v 
and v g , where v g = (l — rj)vK(j>- We initialize the integrations 
at the mid plane with y — z — v z = 0, and x — 1 A.U. 
The in-plane velocities are initialized according to equations 
9 and 10, i.e. v y = v r and v x = v#. The equations are 
integrated for each particle until it reaches 0.1 AU in semi¬ 
major axis, or until the integrator has taken four million 
steps. 

We carry out several hundred integrations 1 with s m 
1 meter, since in the fiducial disk model these are the par¬ 
ticles that traverse the r 3 — 1 line (see figure 1) from the 
weakly- to the strongly-coupled regime. We carry out three 
different sets of integrations with different initial values of 
^target, shown in the three panels of Figure 6. The orien- 


1 We employ the publicly-available bsint package available from 
https://github.com/alrexrudy/bsint. We have verified that it pre¬ 
serves orbital parameters to ~ 10“ 11 over 3000 orbits for our 
setup without the drag and Magnus forces, i.e. the two-body prob¬ 
lem. 


tation of the spin axis Cj is chosen from the uniform dis¬ 
tribution over the surface of the unit sphere. In Figure 6, 
the color shows the z-component of this orientation, Cj z , so 
that red lines show spins aligned with the orbital angular 
momentum vector, while blue lines show anti-aligned runs, 
and lightly-colored lines have a spin-axis nearly lying in the 
orbital plane. 

In all cases the spin has a noticeable effect on the par¬ 
ticles’ trajectories, conceivably extending the particle’s life¬ 
time in the disk by factors of 2 if the spin orientation re¬ 
mained constant. We see that the inward drift velocity of 
the particles is basically determined by cj z , with spin-aligned 
particles flowing in fastest, and progressively more anti¬ 
aligned particles flowing in more slowly. This is an intuitive 
result, since the Magnus force for aligned particles will point 
towards the star, increasing the circular velocity it needs to 
stay in orbit at a given semi-major axis. This increase in the 
velocity will increase the drag force, and shorten the inflow 
time. 

Figure 7 shows more details of the S = 1 run. From top 
to bottom, the panels are inclination, velocity relative to the 
gas, the ratio of the lift to the drag force, the dimensionless 
stopping time, and the fraction of the relative velocity in 
the azimuthal direction. In this more detailed view of the 
integrations, Cj z once again determines the trajectory of the 
particle. For many of the quantities, this is simply a result 
of the aligned (red) particles reaching the center of the disk 
more quickly, with the anti-aligned (blue) particles behav¬ 
ing quite similarly later on. To account for this common 
structure of the trajectories, we can plot each quantity as a 
function of semi-major axis instead of time, averaging over 
some time period to remove the oscillations visible in Figure 
7. 

The result is shown in Figure 8, where each quantity has 
been averaged over a moving 1-year window. Qualitatively 
the paths of the particles do get closer together. However, 
many of the quantities still display substantial differences 
as a result of the particles’ different spin orientations. At a 
fixed radius, the aligned particles have higher inclinations, 
velocities relative to the gas, and Fl/Fp. 

Perhaps most dramatically, the aligned particles experi¬ 
ence negative values of v^—Vgas- In other words, the particle 
experiences a tailwind, yet still spirals in towards the cen¬ 
tral star. In fact, these particles spiral into the star more 
quickly than an analogous non-spinning particle. We can 
also clearly see a qualitative change in the behavior of the 
particles around a = 0.4, particularly the spin-aligned par¬ 
ticles. At this point, meter-sized particles are crossing in to 
the high-Re regime, wherein Fl/Fp approaches its maxi¬ 
mum value (see equation 18). We will investigate this further 
analytically in the next section. 


6 IMPLICATIONS FOR PLANETESIMAL 
FORMATION 

In the previous section we saw a variety of intriguing phe¬ 
nomena in our numerical simulations. Here we attempt to 
understand these results in more detail and their implica¬ 
tions for planet formation scenarios. Examining Figure 8, 
we see that the spread in velocities, like essentially all of the 
time-averaged quantities, is determined by uj z . Moreover, the 
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Figure 9. The vector field T as a function of v re l in the frame co-rotating with the gas. T is evaluated taking u z = +1, s = 1 meter, 
Pgrain = 1 g cm -3 , and M* = 1 Mq. The left column shows a scenario with strong spin, S = 1, while the right column shows S = 0.01. 
Each row shows a different heliocentric radius. The root of T is marked with a red symbol, while the root of the analogous vector field 
evaluated with Cj z = — 1 is shown in blue. The difference in velocity between these two roots in each panel accounts for most of the 
velocity dispersion seen in the numerical integrations. We also show for reference a red dashed line representing a circle of radius ijvk 
centered at zero relative velocity. The color of the streamlines denotes the timescale for a particle to move a distance t\vk in this space. 
The timescales are typically a few orbital times or less, meaning that particles can quickly reach these equilibrium points. 


spread in averaged relative velocities increases substantially 
as the particles make the transition from the weakly-coupled 
t s > 1 to the strongly-coupled r s < 1 regime, but then de¬ 
creases again at small heliocentric radii. In the following we 
show that these results can be explained quantitatively by 
finding the force-free (equilibrium) velocities in a co-rotating 
reference frame. We then use this result to extend our sim¬ 
ulation results to a much wider variety of particle sizes and 
heliocentric radii. 


6.1 Analysis in the co-rotating frame 


In a rotating frame, the force on a particle can be written 


GM*m p x 
(r 2 + « 2 ) 3/2 


+ Fd + Fl + Act 


(38) 


where x is the position of the particle, and v is its velocity. 
The cylindrical radius is r = y/x 2 + y 2 . The last term repre¬ 
sents the “fictitious” Coriolis, centrifugal, and Euler forces, 


Act — —2 m p Q Xv — m p Vl x (Q x x) — x x (39) 

The drag and Magnus forces are of course dependent on 
the velocity and other variables. If we narrow our focus to a 


frame co-rotating with the gas, i.e. Q = (1—tj)Qkz, then v = 
v re l . This is convenient because both of the hydrodynamic 
forces depend directly on v re i . For simplicity we also restrict 
ourselves to z = 0 and co x — uj v — 0 , i.e. the spin of the 
particle must be perfectly aligned or anti-aligned with the 
orbital angular momentum. We know from the results of 
our numerical simulations that for moderate spins, these two 
extremes more or less bracket the particles with intermediate 
alignments. 

In this 2D approximation, we can write out each com¬ 
ponent of the force, 


Tr = 


T& — 


GM*m p v r 

—77(2 — 77 ) --- F d - 

T Z Vrel 

—Fl —-—b 2m p Qic(l — 

Vrel 

-Fd-^~ + Fl— -2ra p Qic(l - f])v r 

Vrel Vrel 

3 / 13 

+ -m p il K v r II - — T) 


(40) 


(41) 


In this context, v re i = yA + v^. The factor of 2 in the 
first term of equation 40 arises from the fact that we are 
in a frame co-rotating with the gas, whose angular velocity 
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Figure 10. The velocity dispersion of particles as a function of radius, normalized to ijvk , the difference between the gas and Keplerian 
orbital velocities (left panel), and the resulting value of the Toomre Q parameter for the dust disk (right panel). Lines are shown 
for different values of the dimensionless spin S. The black, red, and blue lines use Equation 42 to estimate a in both panels, i.e. the 
velocity dispersion is assumed to arise purely from the difference between the equilibrium v re \ of particles aligned or anti-aligned with 
the disk angular momentum. The light blue line with error bars shows the velocity dispersion of the ensemble of 1-meter particles whose 
trajectories were numerically integrated in Section 5. This velocity dispersion is close enough to the approximation of Equation 42, that 
we will adopt a a (which is easy to calculate) throughout the rest of this work. 


is a factor of 1 — r/ slower than Keplerian. This factor gets 
squared when evaluating the centrifugal force. The factor of 
13/21 in the final term comes from the dQ/dt term in the 
Euler force, and the resulting radial derivatives of Qk and 
r] following an application of the chain rule. 

Figure 9 shows streamlines of this vector field as a func¬ 
tion of v r and v# for various locations in the laminar disk 
model (different rows) for a large and a small value of the 
dimensionless spin S — suj/v re i (left and right columns). 
The vector field for Cj z > 0 is shown, while the field for 
<jj z < 0 is not. Each panel of Figure 9 has a red marker de¬ 
noting T r — Tcf — 0. The red dashed line indicates a circle 
of radius t)Vk centered on v r — — 0. In the standard 

no-spin case, we would expect roots of T to lie near this line 
in the weakly-coupled regime, and as a particle entered the 
strongly-coupled regime the root would move towards zero 
velocity, which we do in fact see. 

In addition to the red symbol, each panel of Figure 9 
includes a blue symbol, which is a root of the T vector field 
(not shown) when the spin is anti-aligned with the orbital 
angular momentum. In the low-spin case, we see that the 
red and blue points are virtually indistinguishable, i.e. the 
direction of the spin (and the spin itself) does not matter a 
great deal. As we expected from the numerical integrations, 
we see that in the high-spin case the zero-force velocities are 
appreciably different. 

We posit that a population of spinning particles would 
have substantial velocity dispersions arising from two dif¬ 


ferent effects visible in our numerical simulations. The dif¬ 
ference in position between the red and blue points in each 
panel of Figure 9 demonstrates that even if the particles are 
nearly in equilibrium (i.e. T r — — 0), particles that are 

identical in every respect except their spin orientation will 
have different velocities. This explains the spread in time- 
averaged values of v re i at fixed radius visible in Figure 8. 

The second effect arises from the fact that T r / T<f> ^ 0. 
This is quite clear in the numerical simulations, given the 
short-period oscillations in properties of the particles’ tra¬ 
jectories (Figure 7). These are the result of the continuous 
oscillations in Cl oc uxv re i that occur as particles with fixed 
spin-axis Co (in the inertial frame) rotate about the star, and 
hence change v re i • This effect tends to be much smaller than 
the spread in velocities owing to the difference in equilib¬ 
rium velocities between spin-aligned and anti-aligned parti¬ 
cles. The exception is for particles whose spin axis lies nearly 
in the orbital plane, in which case Cj x v re i has order unity 
oscillations. 

The streamlines in Figure 9 are colored by a quantity 
similar in spirit to the dimensionless stopping time, namely 
m p rjVK^K/T This is the ratio of the time it would take 
a particle to traverse a distance t\vk in this space to the 
orbital time. This value tends to be of order unity or less, 
meaning if the particle is out of equilibrium it will move a 
substantial distance towards equilibrium in a single orbit. 
Moreover the streamlines lead more or less directly to the 
equilibrium point for the small-radii cases, indicating that 
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the particle should be found close to its equilibrium point. If 
the Magnus force were stronger relative to the drag force, the 
streamlines would circle the equilibrium point many times, 
meaning that the particle could plausibly be frequently out 
of equilibrium, however this scenario seems unlikely based 
on equation 18. 

Since we expect the particles to be near equilibrium, i.e. 
T r — T<$, — 0, we have a reliable way of estimating v r and 
Vcf) for Co — ±z. Namely, we can numerically find the root of 
F(v) for arbitrary values of the particle size, location in the 
disk, and dimensionless spin. We denote these two equilib¬ 
rium velocities u+, and V-. We posit that the difference in 
these two velocities is a reasonable proxy for the velocity dis¬ 
persion of particles in the disk, to within a factor of a few. 
We find that the results of the numerical integrations are 
well-matched by the following crude estimate of the velocity 
dispersion 

a = ^\v+-V-\. (42) 

Figure 10 shows a/r/VK as a function of radius for several 
different spin values (left panel). When S = 1, we can di¬ 
rectly compare this prediction to the standard deviation of 
v re i for our randomly-oriented sample of particles from the 
numerical simulations in the previous section. This a nor¬ 
malized by rjVK is shown as the light blue points with error 
bars, and it indeed agrees reasonably well with our estimate 

а. Since it is not numerically difficult to calculate v ±, equa¬ 
tion 42 allows us to estimate the velocity dispersion for any 
particle size, spin rate, and heliocentric radius. It would also 
be straightforward to modify any ingredient of the fiducial 
model, e.g. the mass of the star, the gas density, or the mean 
density of the particles. 

б. 2 Gravitational and Streaming Instabilities 

The velocity dispersion of particles plays an important role 
in planetesimal formation. In a scenario like that envisioned 
by Goldreich & Ward (1973), the stability of the dust disk 
is determined by the Toomre Q parameter, which is directly 
proportional to the velocity dispersion at a given heliocen¬ 
tric radius and particle surface density. Q = Qa / (7rGE s ) is 
shown in the right panel of 10. Another common metric is 
the ratio of the vertical component of the star’s gravity to 
the self-gravity of the particle disk (Garaud &; Lin 2004). 
This stability criterion is 


^disk 4jrG(pd ust + Pg ) 

Note that pdust is the mass density of solid material in a 
large fixed region, as opposed to p S: the average density 
of a single grain. Once again using a as a proxy for the 
velocity dispersion of solid grains, we can estimate pdust ~ 
while p g is taken to be the midplane density given 
in equation 1. 

The value of this stability parameter is mapped out in 
Figure 11 as a function of r and s at fixed S , and as a func¬ 
tion of S and s at fixed r. Throughout both panels we have 
kept p s , o — F — Z re \ — 1 and M* — 1M©. As in previous 
figures, this ratio quickly drops to zero at sufficiently large 
radii, where particles are subject to Epstein drag and we pre¬ 
sume the Magnus force to be inoperative. Nonetheless, we 


see that for small enough heliocentric radii and large enough 
dimensionless spin rates, the particles have a sufficient ve¬ 
locity dispersion to keep g*/gdisk > 1, preventing runaway 
vertical gravitational collapse. 

For our default choice of parameters, the critical size for 
a collection of spinning particles to avoid runaway collapse 
is of order s = 10 cm. This is the same scale on which grain 
collisions cease to be effective at increasing the mass of par¬ 
ticles (Blum & Wurm 2008). One could therefore imagine a 
scenario in which particles coagulate up to 10 cm size-scales, 
and these particles collapse via gravitational instability to 
form km-size planetesimals only in regions where the Mag¬ 
nus force is irrelevant. This would pose a problem because 
terrestrial planets are incredibly common at small heliocen¬ 
tric radii (e.g. Petigura et al. 2013). 

This problem has several plausible solutions. For in¬ 
stance, planetesimals might only form at large radii, with 
the resulting planets migrating through the disk to their 
presently-observed locations at small radii. Another possi¬ 
bility is that the collective gravitational collapse may occur 
when particles are appreciably smaller than 10 cm, in which 
case the Magnus force is unlikely to be strong enough to 
prevent this collapse. This is plausible as long as the stick¬ 
ing efficiency of particles is small enough that by the time 
particles can settle to the mid plane, they have not already 
coagulated to sizes larger than ~ 10 cm (Goldreich & Ward 
1973; Garaud &; Lin 2004). It is also possible that the typical 
spin of particles is simply not large enough for the Magnus 
force to be important at small radii. 

While one or more of these possibilities is likely true, an 
even easier solution is to pass the buck. The laminar dust 
disk assumed in the gravitational collapse scenario is eas¬ 
ily disrupted by even moderate turbulence in the gas disk, 
and even in laminar disks, the dust disk is subject to Kelvin 
Helmholz instabilities from vertical shear (Weidenschilling 
1980). The most popular mechanism to concentrate the 
dust enough for it to collapse gravitationally is the stream¬ 
ing instability. Youdin & Goodman (2005, YG05 hereafter) 
showed that the two-fluid gas+dust equations of motion are 
unstable to axisymmetric perturbations, even in the absence 
of self-gravity and vertical stratification. If the particles are 
allowed to move with respect to the fluid, i.e. r s > 0, these 
instabilities grow and concentrate particles on timescales 
slower than the dynamical time, but quickly enough that 
the particles do not drift into the star. The ability of the 
Magnus force to affect planet formation therefore likely de¬ 
pends on its effect on the streaming instability. 

We find that the Magnus force may indeed stabilize 
certain modes of this instability, though the region of pa¬ 
rameter space where this is true may be fairly small. If the 
orientation of particle spins is approximately isotropic, we 
can estimate that the diffusion coefficient associated with 
the Magnus force is 

D ~ crL, (44) 

where once again we take a « a. L is the effective mean free 
path of the particles. To estimate this, we need to make a 
strong assumption about how frequently a particle’s spin is 
reoriented. If collisions are responsible for keeping the parti¬ 
cles spinning, we expect that the mean free path will be the 
particle’s velocity times the typical time between collisions. 
We also expect that L < H, the scale height of the particle 
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Figure 11. Slices of the ratio of g*/gdisk 5 a criterion used for the stability of a dust disk to gravitational collapse. Values greater than 
1 mean that the disk is stable to collapse. White contours show g*/gdisk = 0.1, 1, and 10. While the red and blue lines show where in 
this parameter space r s = 1, depending on whether their spins are aligned or anti-aligned with the gas angular momentum. This figure 
assumes that the velocity dispersion of the dust is well-approximated by Equation 42. To evaluate this equation, we also assume that 
p s = 1 g cm -3 , F = Z re i = 1, and M* = 1 Mq. If monolithic gravitational collapse is relevant in disks, this figure shows that the Magnus 
force alone is strong enough to stop it for large enough dimensionless spin S and particle size s, and small enough heliocentric radius r. 


disk. We therefore take 


L — min ( at c , H ), 


(45) 


with t c given by equation 33. 

Following section 3.2.2 of YG05, we adopt H ~ rjr, and 
we estimate the diffusion time for modes of wavenumber k , 


tMD — 


4tt 2 

k^D' 


(46) 


The subscript MD is meant to distinguish diffusion owing to 
the Magnus force from other sources, e.g. turbulence. The 
largest modes we expect to be physically relevant will be of 
order the particle scaleheight, in which case k ~ 2n/H , but 
substantially smaller modes tend to grow quickest. YG05 
provide the following fit 2 to the transverse wavenumber that 
maximizes the growth rate of the streaming instability at 
fixed r/rkz — 1 


r\rk x - 


f (2 t s / 9 3 )- 1/2 

\vWsU oa 


fg > 1/2 
fg < 1/2- 


(47) 


Here f g = p g /(p g + Pdust) is the gas fraction. For the r s <C 
1 limit in which the derivation is valid, the dimensionless 
wavenumbers prk x can range from order unity to hundreds. 
The corresponding growth rate of the streaming insta- 


2 The conditionals f g > 1/2 and f g < 1/2 appear to be erro¬ 
neously reversed in Equation 32 of YG05. 


bility is approximately (equation 44 of YG05) 

1 /tsi = Mpfgifp - fa) 2 t s Q k, (48) 

where f p = 1 — f g is the particle density fraction of the 2-fluid 
system. Using prk z = 1 and k x as defined by equation 47, we 
can map the ratio of the streaming instability growth rate to 
the diffusion rate, tMD It si- The diffusion rate is evaluated 
at a wavenumber larger than k = \Jk% + k\ by a factor of 2i r 
to correct for the fact that the wavenumber corresponding to 
one scale height is 27r/(?7r), not l/(rjr). Figure 12 shows this 
ratio for a pair of slices through parameter space. We caution 
that this ratio is only an order of magnitude estimate, since 
it is unclear precisely which wavenumbers to use, how to 
define the diffusion coefficient, and since the growth rate 
from equation 48 is derived from an approximate dispersion 
relation which is not applicable for r s near or greater than 
unity. Moreover there is some ambiguity about which r s to 
use, since it varies as a function of S. 

The only region of the parameter space where we find 
that the Magnus force acts more quickly to diffuse parti¬ 
cles than the streaming instability acts to aggregate them 
occurs near r s — 1, precisely where Equation 48 becomes 
unreliable. This is not a coincidence, as both processes have 
their largest effect near r s — 1. The fact that diffusion ow¬ 
ing to the Magnus force is a strong enough effect that it 
is comparable in magnitude to the streaming instability in 
this critical regime suggests that further investigation in the 
non-linear regime is warranted. 

Over the past 10 years, numerous numerical sim- 
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Figure 12. Slices through parameter space of the streaming instability growth timescale tsi compared to the diffusion timescale owing 
to the Magnus force tMD • White contours are shown at tsi/^MD — 0.1, 1, and 10. Contours where r s = 1 are shown in red and blue for 
particles with spins aligned and anti-aligned with the gas angular momentum respectively. In both panels S = 1, and other parameters 
of the problem are kept fixed at their fiducial values (e.g. F = Z re i = 1). There are parts of the parameter space where tsi/^MD >1, 
meaning that the fastest growing mode of the streaming instability may be suppressed by diffusion owing to the Magnus force. 


ulations have been performed to study the nonlinear 
growth of the streaming instability (e.g. Youdin h Johansen 
2007; Johansen & Youdin 2007; Johansen et al. 2007; 
Yang & Johansen 2014). Even if sufficiently small or suffi¬ 
ciently slowly spinning dust grains are initially subject to the 
streaming instability, at some stage during the collapse and 
coagulation of the concentrated dust the Magnus force may 
become important. As the particles increase in size, their 
stopping times and Reynolds numbers may pass through 
the critical values of r s ~ 1 and Re ~ 800. Indeed, within 
r < 1 A.U. particles pass through these critical values simul¬ 
taneously when the particles are roughly 1 meter in radius, 
leading to velocity dispersions of order t]Vk- Moreover, the 
spin of particles in these clumps may be large, since collisions 
occur frequently. Whatever small rotation rate is present in 
the initial large-scale region will be amplified by the conser¬ 
vation of angular momentum during the collapse to smaller 
size scales. Nesvorny et al. (2010) showed that this mecha¬ 
nism could account for binaries in the Kuiper Belt, and we 
suggest that it could also lead to high spin rates for individ¬ 
ual objects formed in this process. 


7 SUMMARY 

We have introduced the Magnus force as a potentially im¬ 
portant piece of microphysics for solid bodies in gaseous pro¬ 
toplanetary disks. Rotating particles, by an asymmetry in 
the wake they leave in the background fluid, experience a 
force roughly perpendicular to both their direction of rela¬ 
tive motion and their spin axis. Just like the drag force, the 


Magnus force depends on the dimensionless parameters of 
the flow, namely the Reynolds number Re, the dimension¬ 
less spin S = suj/v re i , and the ratio of the gas mean free 
path to the particle size, A/s. This dependence produces 
three (particle-size-dependent) regions in the disk. 

In the outermost region, we expect that the Magnus 
force will be irrelevant because A >s. As particles drift in¬ 
ward in the disk eventually the gas begins to behave as a 
fluid, but the Magnus force remains comparatively weak. 
The only exception is if S 1, which we consider improb¬ 
able. In the inner disk when Re >800, the Magnus force 
becomes comparable to drag if S ~ 1. 

Clearly a great deal depends on the spin rate. In the 
modern solar system, most of the available data on spin 
rates is for kilometer-sized objects likely held together by 
self-gravity and confined to moderate spin rates. The small¬ 
est objects for which data is available show a hint towards 
higher spin-rates, scaling roughly as u oc s~ , in which case 
S oc v~^. In protoplanetary disks, spins can be both induced 
and damped by hydrodynamic drag, leading to a median 
S > 10 -3 with an extremely broad log-normal distribution. 
Meanwhile frequent collisions between particles of compara¬ 
ble size may keep a large population of particles spinning 
quickly near S ~ 1. Notably the orientation of the parti¬ 
cle spins in these scenarios is likely close to isotropic. The 
coagulation of large groups of particles via the streaming 
instability or gravitational collapse may also induce large 
spins during this critical phase where particles may traverse 
the meter barrier. 

Presuming that at least modest spin rates can be sus- 
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tained in the disk, we carried out direct numerical integra¬ 
tions of the equations of motion for individual meter-sized 
objects to understand their dynamics when the Magnus 
force is included. Even when the Magnus force is comparable 
to the drag force, and even when the orientation of the parti¬ 
cle is assumed to remain constant throughout the particle’s 
life, the effect on particle lifetimes in the disk is moderate, at 
most a factor of two. Qualitatively, the spinning particles be¬ 
have quite similarly to the non-spinning particles, gradually 
spiraling in to the central star. The trajectory followed by 
the spinning particles is determined by the component 
of the spin rate parallel to the gas angular momentum, with 
aligned particles drifting in faster, and anti-aligned particles 
surviving longer relative to their non-spinning counterparts. 

The Magnus force can induce small but non-negligible 
out-of-plane motions, and perhaps most surprisingly, it can 
reverse the azimuthal velocity of particles relative to the 
gas. Non-spinning particles always orbit more quickly than 
the pressure-supported gas, but particles whose spins are 
aligned with the gas angular momentum actually orbit more 
slowly than the gas. This can be understood geometrically 
in a frame co-rot at ing with the gas. When a particle enters 
the strongly-coupled (r 8 < 1) regime, its inward radial mo¬ 
tion becomes larger relative to its azimuthal velocity, and so 
the drag force begins to point outward. Spinning particles 
of opposite alignment will have Magnus force vectors point¬ 
ing ±90° away from the drag force vector. Particles with 
the Magnus force pointing in the negative azimuthal direc¬ 
tion will then acquire negative equilibrium velocities in this 
frame, i.e. they will orbit the star more slowly than the gas! 

At a fixed radius, the velocity of particles with different 
spins can vary by up to a factor of two. Although a fac¬ 
tor of two is unlikely to be significant in terms of particle 
lifetimes in the disk (the meter barrier is a problem regard¬ 
less of whether particles spiral inwards in 100 or 200 years), 
this spread in velocities is large enough to have implications 
for theories of planetesimal formation. We show, by analyz¬ 
ing the equilibrium velocities of particles with oppositely- 
aligned spins in a frame co-rotating with the gas, that the ve¬ 
locity dispersion of a population of spinning particles is large 
enough to prevent monolithic runaway gravitational collapse 
of the dust layer in the inner regions of the disk for parti¬ 
cles larger than about 10 centimeters. This suggests that 
the planets found in abundance at small orbital periods in 
observations either formed at larger radii, or that the plan- 
etesimals which would go on to form these planets formed 
from dust grains smaller than 10 centimeters. If this is the 
case, the gravitational collapse that formed these planetesi- 
mals had to occur more quickly than the timescale for grains 
to coagulate to sizes larger than 10 centimeters. In other 
words, the particles are subject to a race between settling 
by vertical gas drag and growth by coagulation. Counterin¬ 
tuitively, if coagulation wins planetesimal formation may be 
suppressed by the diffusive effects of the Magnus force. 

Monolithic gravitational collapse can also be disrupted 
by larger-scale turbulence in the disk and shearing instabil¬ 
ities that tend to stir up very thin dust disks. Attention has 
therefore turned to other means of concentrating particles, 
especially the streaming instability. The two-fluid gas plus 
dust equations of motion with non-zero r s exhibit a grow¬ 
ing axisymmetric mode which concentrates dust on reason¬ 
ably short timescales. We compare the linear growth rate 


with the timescale on which particles diffuse owing to the 
Magnus force, assuming the spin axes of the particles are 
reoriented on a collisional timescale. We find that diffusion 
owing to the Magnus force is comparable in strength to con¬ 
centration owing to the streaming instability when r s ~ 1, 
although this result will require multi-dimensional simula¬ 
tions to verify. The Magnus force may play an even more 
important role in the nonlinear evolution of this instability 
since the collapse of large-scale modes will tend to spin up 
particles by conservation of angular momentum. 
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